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ABSTRACT 

Large scale magnetic field threading an accretion disk is a key ingredient in the jet formation model. The 
most attractive scenario for the origin of such a large scale field is the advection of the field by the gas in the 
accretion disk from the interstellar medium or a companion star. However, it is realized that outward diffusion 
of the accreted field is fast compared to the inward accretion velocity in a geometrically thin accretion disk if the 
value of the Prandtl number P m is around unity. In this work, we revisit this problem considering the angular 
momentum of the disk is removed predominantly by the magnetically driven outflows. The radial velocity 
of the disk is significantly increased due to the presence of the outflows. Using a simplified model for the 
vertical disk structure, we find that even moderately weak fields can cause sufficient angular momentum loss 
via a magnetic wind to balance outward diffusion. There are two equilibrium points, one at low field strengths 
corresponding to a plasma-beta at the midplane of order several hundred, and one for strong accreted fields, 
/3 ~ 1. We surmise that the first is relevant for the accretion of weak, possibly external, fields through the outer 
parts of the disk, while the latter one could explain the tendency, observed in full 3D numerical simulations, of 
strong flux bundles at the centers of disk to stay confined in spite of strong MRI turbulence surrounding them. 
Subject headings: accretion, accretion disks, galaxies: jets, magnetic fields 



1. INTRODUCTION 

Jets/outflows are observed in different types of the sources, 
such as, active galactic nuclei (AGNs), X-ray binaries, and 
young stellar objects, which are probably driven from the ac- 
cretion disk through the magnetic field lines threading the 
disk (see reviews in ISpruitl 119961: iKonigl & Pudritzl 120001: 
iPudritz et al.1 12007b ISpruitl 1201 Oh . The large scale magnetic 
field co-rotates with the gases in the disk, and the jets/outflows 
are powered by the gravitation energy released by accretion of 
the gases through the ordered field threading the disk. A large 
scale magnetic field, of uniform polarity threading the inner 
parts of the disk is probably a key ingredient in this jet forma- 
tion model. 

The origin of such a field is not well understood, however, 
since the net magnetic flux threading a disk cannot be pro- 
duced or changed by internal processes in the disk (a conse- 
quence of the solenoidal nature of the magnetic field, see e.g. 
Spruit 2010). A net magnetic flux in the inner disk must there- 
fore either be inherited from initial conditions, or somehow be 
accreted from a larger distance; ultimately for example from 
the interstellar medium or a companion star (cf. Bisnovatyi- 
Kogan & Ruzmaikin 1974, 1976). 

One could imagine a steady state in which the inward 
advection of the field lines is balanced by the outward 
movement of field lines due to magnetic diffusion. In 
conventional isotropic idealizations of a turbulent plasma, 
it is expected (e.g. Parker 1979) that v ~ r\ ~ /v t , in 
which / is the largest eddy size, and v t is turnover ve- 
locity, i.e. P m ~ 1, Whether this is actually the case in 
MRI turbulen ce has been investigated using numerical sim- 
ulations (e.g.. lYousef et al.1 120031: IL esur & Longarettij 120091: 
iFromang & Stand 120091: iGuan & Gammid 120091) . The re- 
sults all suggest that the effective magnetic Prandtl number 
is around unity. Fromang et al. (2009), for example, measure 
P m w 2. For such Prandtl numbers a vertical field (perpendic- 
ular to the disk plane) can indeed be dragged efficiently by an 



advection dominated ac cretion flow dCao 201 1), which is hot 
and geometrically thick ( iNaravan & Villi 9941, 1 1995b . 

The advection of the field in a geometrically thin (H/r <C 
1 ), turbulent accretion disk is inefficient, however, because the 
radial component of the magnetic field diffuses much faster 
across the disk, on a time scale ~ H 2 /rj. As a result, a mag- 
netic field with inclination B r /B z ~ 1 actually diffuses out- 
ward on a time scale of order H/r shorter than a purely verti- 
cal field (van Ballegooijen 1989). An inclined field is a neces- 
sary consequence of effective accretion of the field, however, 
since the accumulation of field lines in the inner disk exerts a 
pressure that causes the field above the disk to spread outward. 
While such an inclined configuration is favorable for laun ch- 
ing a flow dBlandford & Pavndl 19821: ICao & Spruitll 1994b . it 
raises the problem how it can be accreted effectively against 
the action of magnetic diffusion. 

Several alternatives were suggested to r esolve the difficulty 
of field advection in thin accretion disks. Spruit & Uzdensky 
(2005) suggested that a weak large-scale magnetic field 
threads the disk in the form of localized patches in which the 
field is strong enough to cause efficient angular momentum 
through a magnetic wind. General relativistic magnetohydro- 
dynamic (GRMHD) simulations of an accretion torus embed- 
ded in a large-scale magnetic field showed that a central mag- 
netic flux bundle, once formed from a suitable initial condi- 
tion, can survive in spite of MRI turbule nce present in the 
disk surrounding it (Beckwith et al. 2009). The calculations 
bv iGuilet & Qgilvid d2012ab and Guilet & Ogilvie (2012b) 
show that the accretion velocity of the gas in the region away 
from the midplane of the disk can be larger than that at the 
midplane of the disk, which may partially solve the problem 
of too efficient diffusion of the field in thin accretion disk. 

The outward diffusion of the field could be balanced by ac- 
cretion if a process can be found that increases the accretion 
velocity by a factor ~ r/H relative to the rate due to the turbu- 
lence alone. We explore here the conditions under this can be 
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achieved by a magnetic wind generated by the weak magnetic 
field that is to be accreted. 

2. MODEL 
2.1. Model assumptions 

Apart from the accreted field, the disk model we use is a 
standard a-disk model, i.e. with a viscosity v parametrized 
as v = ac s H, where c s is the (isothermal) sound speed, H is 
the scale height of the disk, and a ~ 0.01 -0.1 (the range of 
values measured in MRI simulations). The field to be accreted 
is assumed to be sufficiently weak that it does not suppress 
magnetororational instability. The disk thus contains an MRI- 
generated field as well as a weaker field of uniform polarity 
threading the disk. 

Since the accretion velocity that is to be achieved exceeds 
the viscous rate, we simplify the analysis by assuming the ac- 
cretion flow to be dominated by the angular momentum loss 
from the accreted field, ignoring the viscous contribution from 
the MRI turbulence. The resulting accretion rate then has to 
exceed outward diffusion due to the MRI turbulence. The tur- 
bulence is assumed to produce an effective magnetic diffusiv- 
ity corresponding to a magnetic Prandtl number P m ~ 1 . 

The model needs a prescription for the angular momen- 
tum loss produced by the accreting weak magnetic field. For 
this we use the Weber-Davis model for a magnetically driven 
wind, in the 'cold' approximation (in which the gas pres- 
sure force is neglected). It leads to a simple description in 
terms of the field strength and mass loss rate (Mestel 2012, 
Spruit 1996). The Weber-Davis model strictly applies only 
to the 'split monopole' configuration, in which the poloidal 
field is purely radial. Its properties are found to be a rather 
good approximation for more general poloidal field shapes as 
well (Anderson et al. 2005), which makes it adequate for the 
present purpose. 

The angular momentum loss in this model can be char- 
acterized by a single dimensionless constant, a 'mass load 
parameter>(c.f. Michel 1969, Mestel, 2012): 



jU = x47rfir /B p 



(1) 



where B p = (B^+Bf) 1 / 2 is the poloidal component of the ac- 
creted field at radius ro on the surface of the disk, and SI the 
rotation rate at ro (cylindrical coordinates r,<fi,z)- x is a con- 
stant along the flow, it is a measure of the 'mass flux per field 
line': 

X = pv p /B p , (2) 

where p, v p , B p are the mass density, poloidal velocity and 
poloidal field strength. It is related to the mass flux per unit 
surface area from the disk »i w by 



= B zX . 



(3) 



The asymptotic velocity of the wind decreases monotonically 
with increasing //, and for // = 1 equals the rotation velocity at 
ro. The angular momentum loss increases monotonically with 
H. 

It turns out that the conditions for effective accretion of the 
field can be satisfied when /i > aP m r/H. This is described 
in the following sections, using a specific model for the disk 
structure. 

2.2. The magnetic field of the disk 



The angular momentum equation for a steady accretion disk 
with outflows is 
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(4) 



where £ is the total surface mass density of the disk (count- 
ing both sides), 7j n the total torque per unit of surface area of 
the disk (counting both sides) due to the magnetically driven 
outflows. Integrating Eq. |@J, yields 



27rrEv r r SI = 27rr^Er 



dr 



-2nf m (r) + C, 



(5) 



where the value of the integral constant C can be determined 
with a boundary condition on the accreting object, and 



rT m 



df m (r) 
dr 



(6) 



The magnetic torque in unit surface disk area exerted by the 
outflows is 



T m = 2ra w r 4 Sl, 



(7) 



where ra w is the mass loss rate in the outflow from a unit 
surface area of the disk (single sided), r A is the (cylindrical) 
Alfven radius of the outflow, Sl(r) is the angular velocity of 
the disk, and r is the radius of the field line foot point at the 
disk surface. The dimensionless mass load parameter fi of the 
outflow is 



47rp w v w Slr Andr 



B P B, 



(8) 



where p H 



v B p /B : is the mass flux parallel to the field 



line. In the cold Weber-Davis model the Alfven radius is 



^(l + M- 2/3 ) 



"I 1 / 2 



(9) 



[F or more detailed discussion of this model see Mestel (2012) 
or Spruit (1996)]. MHD simulations of axisymmetric mag- 
netically driven flows have shown that relations like © and 
( [Tol l below are fair approximations for more general config- 
urations than the Web er-Davis 'split monopole' (cf. Fig. 7 in 
lAnderson etaTI 120051) . Substituting Eqs. © and ® into 0, 
we find 
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rB^(l + M - 2 / 3 ). 



(10) 



The radial velocity of an accretion flow in which the angular 
momentum is removed predominantly by the outflows can be 
estimated from Eq. (0): 



v r • 



7m 



n -l 



dr 



(r 2 Sl) 



2T m 



(11) 



where we have assumed that the rotation is approximately Ke- 
plerian, SI w S1k- [This is sufficient for the following esti- 
mates, but has to be made more pre cise when considering the 
wind launching conditions (section l23l ).l The mass accretion 
rate of the accretion disk is 



M = -2irrY,v r ■ 



4tt7™ 
SI 



(12) 



where Eq. dTTb is used. We can compare the mass loss rate in 
the outflows with the accretion rate through the disk. Using 
©, (OB and O: 



d\\\M 4Trr 2 m w 1 
dlnr M 3 



(13) 
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This shows that for large mass loading parameters, // > 1, a 
fraction ~ 1/3 of the accretion rate can be lost in the wind, if 
the wind is sustained across the disk over a distance of order r. 
For n < 1 the mass loss in the wind is not significant compared 
with the accretion rate. Substituting ( [TUt into ([TO . we obtain 



6c s g(fi) 

Vr ~ n H p p Cs 



(14) 



where g = /i(l + /i 2 / 3 ), c s is the sound speed of the gas in the 
midplane of the disk, and 
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(15) 



is a measure of the poloidal field strength at the surface rela- 
tive to p c , the gas pressure at the midplane of the disk. The 
magnetic field is advected inwards on a timescale 



7"adv 



r (lHP p 



(16) 



In the limit H/r -C 1, the radial component of the field dom- 
inates the outward diffusion timescale t^s of the magnetic 
field (van Ballegooijen 1989, Lubow et al. 1994). Let 



Ko=B z /B r , s 



(17) 



be the inclination, with respect to the horizontal, of the ac- 
creted field at the disk surface (s). The diffusion time scale is 
then of order 

Tdif~ , (18) 

V 

where rj is the magnetic diffusivity. [This dependence holds 
as long as 1/kq > H/r; in the opposite case of a nearly 
vertical field B r /Bo < H/r, the diffusion time is of order 
r 2 /i]]. The magnetic Prandtl number is defined as P m = -q/v, 
where v is the turbulent viscosity. With the conventional a- 
parametrization, v = ac s H we can then write Eq. ( fT8l as 



Tdif 



rn 



(19) 



For a steady state, the advection of the field in the disk has 
to balance the diffusion, i.e. r a d v = Tdif. With ( fToT ). ( fT9l this 
yields a condition on the mass flux parameter fi: 



Md + M" 2/3 ) : 



aQHfipPn 
6c s k 



which reduces to 



-2/3) _ Q A> P n 

6k 



(20) 



(21) 



in weak field approximation, c s = VLH. If a significant mass 
flux, n ~ 1 can be launched, and assuming a = 0.01, P m = 1, 
and a field inclination of 30° to the vertical, this shows that 
a weak field with j3 ~ 10 3 can still be accreted through the 
angular momentum loss of the wind associated with it. In the 
following sections we investigate this with specific models for 
the structure of the disk. 

2.3. The mass flux in the wind 

The mass loss rate in the outflow is governed by the gas den- 
sity at the position of the sonic point (in the following labeled 
with index s ), where the flow speed equals the sound speed 



(more accurately, the slow mode cusp speed). The mass flux, 
parallel to the field, is thus approximately pc s . s , and per unit 
of surface area parallel to the disk surface, the mass flux is 



B z 
B p 



(22) 



where B p = {B 2 +B 2 ) l l 2 as before. As mentioned, we consider 
only the wind associated with the accreted field. The MRI 
turbulence also produces a field, but being a local process, we 
assume that it does not produce a large scale field that would 
be significant for wind-driven angular momentum loss. 

The mass flux depends sensitively on the surface temper- 
ature T s . Assume a radiative disk, i.e. the vertical energy 
transport through the disk is by radiation. In the diffusion 
approximation, the surface temperature is then related to the 
temperature T c of the disk at the midplane of the disk by 



3t 



(23) 



where r is the optical depth of the disk in the vertical direc- 



tion. In terms of the sound speeds c s 
the midplane, 

4 \ V8 
3r, 



c s c at the surface and 



(24) 



Along a field line that corotates with its footpoint in the disk 
the outflow is governed by an effective potential "J^ff, 



* e ff(r,z): 



GM 



(r 2 +z 2 )!/2 2 



-tfr 2 , 



(25) 



where ft is the angular velocity of the footpoint at the disk 
surface. In this expression and in the following, the (small) 
difference between the rotation rate of the field line and the 
Keplerian value £1k has to be included consistently, since it 
has a strong effect on the launching conditions for the wind, 
through its effect on \t e ff (Ogilvie and Livio, 2001). This dif- 
ference, due to the magnetic stress exerted at the surfaces fol- 
lows from the radial equation of motion, 



rnl-rfl 2 



B 1 



/>'■>/>' 



(26) 



where S the total (two-sided) surface density S w 2p c H. This 
yields 



n = n K 



1- 



2r 



(3 z k H r 2 n\ 
which can be approximated as 



1/2 



n = n K 



l- 



H 2 

r (3 z k 



1/2 



(27) 



(28) 



in weak field case, where (3 Z = %np c /B 2 = /3 p (l + Kq)/kq is the 
plasma-beta of the accreted field evaluated at the midplane ( c ) 
of the disk. 

Next, we make the model more specific and simplify it a 
bit by approximating the sound speed in the launching region 
as constant with height. In addition we take the midplane 
temperature as representative for the interior of the disk. This 
is sufficient for the evaluation of quantities like the rotation 
rate correction in (128V 

The location of the sonic point is close to the maximum of 
the effective potential. The mass flux is determined by the 
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density at the sonic point; a fair approximation for this is to 
treat the subsonic region as if it were in hydrostatic equilib- 
rium. This yields the following estimate for the mass loss rate 
in an isothermal outflow (per unit of disk surface area), 



«0 



(l + K 2)l/2 



A,c s ,s exp [-(* ef f,s - *eff,0)/Cs, s ] , (29) 



where po is the density of the gas in the base of the outflow 
(still to be specified), and ^eff.s and V&eff.o are the effective 
potential at the sonic point and the footpoint at the disk sur- 
face, respectively. The factor involving kq is equal to the ratio 
B z /B ns in Eq. Q/7}. 

The vertical structure of an isothermal disk can be calcu- 
lated with the vertical momentum equation, 



ci^— =-p(z)Sliz- 



dz 



B r (z) dB r (z) 
4ir dz 



(30) 



In principle, the field line shape is computable by solving the 
radial and vertical momentum equations with suitable bound- 
ary conditions (for a detailed discussion see Cao & Spruit 
2002, hereafter CS02). For the isothermal case, an approx- 
imate analytical expression is proposed for the shape of the 
field lines in the flow: 



H 



-vf+vh 2 H~ 



-2)1/2. 



H 



(1 



■ri) 



1/2 



(31) 



where r; is the radius of the field line footpoint at the mid- 
plane of the disk, and 77; = tanh(l) (see CS02). This expres- 
sion reproduces the basic features of the Kippenhahn-Schluter 
model for a sheet of gas suspended again st gravity by a mag- 
netic field ( Kippenhahn & S chliiteJ 1 95 7h . for weak as well as 
strong field cases. As done in CS02, we use a fitting formula 
to calculate the scale height of the disk in the rest of this work, 



H 



Ac 2 



where 



Solving for H: 



H = 



1/2 



B 2 



1 



2(1 -e- l / 2 ) KQ 4irprHn 2 K K ' 



Ok" 



(l-e- 1 /2) / 3 p (l + K 2 ) 



1/2 



which requires 



or 



Pp ft),min — 



ft > ft.n 



(l-e- 1 /2)(l + K 2 ) ' 



(1- 



(32) 



(33) 



(34) 



(35) 



(36) 



The square bracket in (I34t gives the magnetic correction to 
the standard relation between disk thickness and the sound 
speed at the midplane. It reduces to unity for a weak field, or 
when the radial component is small compared with the verti- 
cal component. For a typical value of kq ~ 1, we find ft must 
be > 1.3. This means that the magnetic pressure can be larger 
than the gas pressure in the disk only if the inclination of the 
field line at the disk surface kq is larger than unity. 



Equation (1321 should be a good approximation for the 
present investigation, especially in the thin outer region of the 
disk which is probably the most critical region for the accre- 
tion of a net magnetic flux. Substituting (l29b into ([SJ, the mass 
flow parameter is 



p = — ^5 — eX P r(*eff,s - 

P 



*eff,0)/Cs, s 



(37) 



To estimate the density po at the base of the flow, we note 
that Eq. d37l ) is applicable only at heights in the atmosphere 
where the field is strong enough, relative to the plasma, to 
enforce corotation so the effective potential is relevant for 
the launching process of the wind. Inside the disk, where /3 > 
1, this is not the case. As base of the flow, where the pressure 
is po, we assume the height where the plasma-beta is of order 
unity, i.e. 

ft EE 87T/7o/flp « 1, (38) 



which determines po if ft is given. For most of calculations 
reported below, ft = 1 is used, some with a lower value 0. 1 . 
The mass flux is then related to the poloidal field strength at 
the disk surface by 



Bi 



P0C s ,s = 



87TC S 



Using ( f24b Eq. d37| i can then be written as 
fir 



(3T/4) 1 / 8 exp [-(* eff , s - *eff,0)/cy 



(39) 



(40) 



The factor in front of the exponential tends to be a large num- 
ber, the exponential itself a small one. To evaluate the ef- 
fective potential as a function of height, field line inclination, 
and the (slightly non-Keplerian) rotation rate ft, the model of 
CS02 is used. It is also used for the optical depth connecting 
the surface temperature to the disk midplane temperature. 

3. RESULTS 

The disk is compressed in the vertical direction by the 
curved magnetic field line, which sets an upper limit on the 
magnetic field strength. We plot the corresponding minimal 
ft as a function of the field inclination kq at the disk sur- 
face in Fig.Q] The curved field line also exerts a radial force 
on the disk against the gravity of the central object, which 
provides an additional constraint on the field strength (Fig. [T] 
show the result for different values of the disk temperature). 
The angular velocity fl deviates from the Keplerian value due 
to the radial magnetic force. Fig. [2] shows the dimensionless 
angular velocity fl/fl^ of the disk as a function of kq and 
e = c 2 J{r 2 fll). 

For given values of the disk parameters, i.e., the viscosity a, 
the temperature parameter in the disk, and the optical depth 
t, the dependence of mass loading p on ft can be calculated 
with Eqs. d20T ) and ( f4Qb respectively. These relations are plot- 
ted in Figs. [3][5]for disk-outflow systems with different val- 
ues of the parameters. In all calculations, P m = 1 is adopted. 
It is found that two branches of solutions usually exist for 
most cases. We plot the solutions of the disk-outflow systems 
in Figs. [6] and [7] The solutions with different values of ft 
(ft = 8irpo/B 2 ) at the base of the outflow for the isothermal 
disk are compared in Figs. |8]and|9] The two branches of so- 
lutions correspond to low mass (// <C 1) loaded outflows with 
strong field strength (low-ft), or high mass loaded outflows 
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with relative weak field strength. The optical depth needed to 
connect surface and internal temperature has been computed 
by including a Rosseland mean opacity and electron scatter- 
ing; it is shown in Fig.fTUl 

For comparison, we have repeated the calculations for a 
uniformly isothermal disk, that is, the temperatures of disk 
and wind are assumed to be the same (Figs.|8]and|9). Though 
this is not very realistic, it gives an impression of the sensitiv- 
ity of the results to the model assumptions made. 




o 



FIG. 1 . — Constraints on the poloidal magnetic field strength in terms of /3- 
(the plasma-beta of the accreted field measured at the midplane of the disk, 
see text). Black line: fi 7 . m j n a s constrained by the vertical pressure exerted by 
the curved field line (Eq. l36t . Color s: th e minimal values of /3 p constrained 
by the radial magnetic force (see Eg. 1271 for different disk temperatures, = 
0.01 (red), 2.5 X 10" 4 (green), and 10" 4 (blue). 
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FIG. 2. — Deviation of the angular velocity of the di sk Q relative to Kep- 
lerian due to magnetic stress, as a function of j3 z (Eq. l27t . Field line incli- 
nation kq= 1 (red), \/3 (green), and 2.5 (blue). Disk temperature parameter 
= 0.01 (solid), 2.5 X 10" 4 (dashed), and 10"* (dash-dotted). 

3.1. Discussion 

The magnetic field is dragged inwards by the accretion disk, 
an outflow is launched by this field, and the radial velocity of 
the accretion disk is determined by the rate of angular mo- 
mentum carried away by the outflows. This loss rate in the 
outflows can be estimated by exploring the launching pro- 
cess of the outflow, which depends sensitively on the field 
strength/configuration and the density/temperature of the gas 



at the disk surface. We have used the cold Weber-Davis model 
for determining the mass load parameter /i of the outflow as a 
function of strength of the accreted field and parameters of the 
disk structure: its optical depth (a radiative disk is assumed), 
a temperature-parameter 6, the a-viscosity, and the magnetic 
Prandtl number P m of the assumed MRI turbulence (cf. Eq. 
l20l i. Balancing the resulting accretion velocity with the out- 
ward diffusion by magnetic turbulence determines the condi- 
tions for existence of a stationary disk-outflow system. Figs. 
[3]|5] illustrate the properties of the solutions. It is found that 
two solution branches exist for all cases. The lower branch 
corresponds to high field strength and low /1, i.e., low mass 
loss rate, the upper one corresponds to low field strength and 
high fi, i.e. high mass loss rate (see Figs. [6] and |7). There 
is an upper limit on the field inclination kq at the disk sur- 
face, which increases with disk temperature. Overcoming the 
deeper effective potential barrier associated with a larger in- 
clination requires a higher internal energy of the gas. The 
maximum inclination thus decreases with increasing opti- 
cal depth t of the disk, as the surface temperature decreases 
with increasing t, (keeping other disk parameters fixed, see 
Eq.|2D. 

We have also explored the sensitivity to model assump- 
tions somewhat with solutions for a more drastically simpli- 
fied case, where the temperature is assumed uniform through- 
out, i.e. the temperature of the outflow is the same as the disk 
temperature. This is shown in Figs. [8] and [9] These also show 
the effect of assuming a lower density of the gas at the base of 
the outflow (/3 S = 0.1). The results are qualitatively similar to 
those with /3 S = 1 . 




10" 1 10° 10 1 10 2 10 3 
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FIG. 3. — Mass loading /iasa fu ncti on of /?-. Solid: as determined from 
the wind launching conditions Eq. 440t Broken: mass loading required for 
effective accretion of the field lines l !20t . for disk viscosity a = 0.1 (dashed)) 
and a = 1 (dash-dotted). The intersection points are possible solutions for the 
stationary wind driven accretion problem. Field line inclinations are Kg = 1.5 
(red), v3 (green), and 2 (blue). The disk temperature parameter = 0.01, 
magnetic Prandtl number P m = 1 . 




FIG. 4.— As Fig.[5]for = 2.5 X 10" 4 . 
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FIG. 5.— As Fig.[3] for = 10~ 
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FIG. 6. — Resulting disk-outflow solutions for different disk temperature 
parameters: = 0.01 (red), 2.5 X 10 -4 (green), and 10^* (blue), and different 
values of the disk optical depth: r = 10 2 (solid), 10 3 (dashed), and 10 4 (dash- 
dotted). The viscosity parameter a = 0.1. 





FIG. 8. — The disk-outflow solutions for vertically isothermal accretion 
disks. The colored lines represent the solutions derived with different disk 
temperature, = 0.01 (red), 2.5 X 10 -4 (green), and 10 -4 (blue), respectively. 
The solutions derived with different ratios of gas pressure to magnetic pres- 
sure at the disk surface are indicated with different line types, S = 1 (solid), 
and 0.1 (dashed), respectively. Viscosity parameter a = 0.1. 




FIG. 10. — The optical depth of the standard thin accretion disks with- 
out magnetic field as functions of radius, where rs = 2GM/c~. The opac- 
ity Ktot = «cs + is adopted, where the Rosseland mean opacity kr = 

5 X 10 24 pr c 1 I~ cm 2 g _1 . The red lines represent the results calculated for 
a black hole with M = IOMq, while the blue lines are for a black hole with 
M = 10 8 Mq. The different line types correspond to the disks with different 
parameters: solid (a = 1 and m = 0.1), dashed (a = 0.1 and m = 0. 1), dash- 
dotted (a = 1 and m = 0.01), and dotted lines (a = 0.1 and m = 0.01). 

The dependence of the solutions on the value of a is shown 
in Fig. [7] The magnetic diffusivity rj is scaled with the tur- 
bulent viscosity v, and therefore the diffusion becomes more 
important for the cases with a higher value of the viscosity 
parameter a. In order to compete with the diffusion of the 
field, a large radial velocity of the disk is required for high a 
cases, which corresponds a high rate of angular momentum 
removal by the outflows. We find that the values of the mass 
load parameter /i of the outflow are systematically higher for 
those derived with a larger a (compare the results in Figs. [6] 
andE). 
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The accretion disk is vertically compressed by the curved 
magnetic field, which sets an upper limit on the field strength 
(see Eqs.[3~5land[36ll. The value of /3 z m ; n only depends on the 
field inclination ko at the disk surface, and /3 ; m i n decreases 
with increasing kq (see Fig. [T). There is a force exerted on 
the disk by the curved field against the gravity of the central 
object in the radial direction, which makes the rotation of the 
gas in the disk be sub-Keplerian. The rotational velocity of 
the disk can be quite low if the field strength is sufficiently 
strong (see Eq.l27land Fig.[2]l. The disk is then magnetically 
supported against gravity. Such configurations are likely to 
be unstable to interchange instabilities, however (Spruit et al. 
1995), which effectively cause the magnetic field to spread 
outward and limit the field strength (as observed in the nu- 
merical simulations of Stehle & Spruit 2001). The require- 
ment that the magnetic force is less than the gravity in the 
radial direction provides an additional constraint on the field 
strength. We find that the constraints almost overlap with that 
constrained by the force in the vertical direction when kq < 3, 
while maximal field strength becomes lower (a larger /3mi n ) 
for a disk with relative high temperature and kq > 3 (see Fig. 

ID- 

From Fig. Q] we see that the magnetic field can be very 
strong, e.g., the magnetic pressure can be more than one or- 
der of magnitude higher than the gas pressure in the disk if the 
field inclination hq is sufficiently large. However, the magnet- 
ically driven outflow will be suppressed if ko is too large, be- 
cause the effective potential barrier becomes extremely deep 
in this case. The results show that /? >0.5 is always satisfied 
in the disk-outflow solutions (see Figs. |6] and |7]i. Note that this 
does not mean the magnetic field cannot be strong, as the disk 
is significantly compressed in the vertical direction, which in- 
creases the density of the disk and then the gas pressure for 
given disk temperature. 

3.2. Stability 

We have calculated only stationary solutions, but their sta- 
bility to time-dependent perturbations can already be guessed 
at by inspection of the intersection points in Fig. [3] Near 
the high-beta (low field) solution the mass loading parameter 
(solid line) is nearly independent of the field strength assumed 
(this is because of the assumed value of the plasma-beta at the 
base of the flow). If the field strength were to decrease (to 
the right of the intersection point), the mass loss would need 
to increase in order to maintain a balance between inward ac- 
cretion and outward diffusion (dashed line). Since the mass 
loading actually does not change much, the angular momen- 
tum loss is insufficient to balance outward diffusion for such a 
perturbation. Outward diffusion will then tend to decrease the 
field strength, providing a positive feedback to the perturba- 
tion. This stationary solution is thus expected to be unstable. 

This mechanism of instability is the same as that identified 
in the linear stability analysis of CS02. The timescale of the 
instability is comparable with the dynamical timescale of the 
disk if the magnetic torque is large, which becomes signif- 
icantly small when the magnetic torque is weak (see CS02 
for the detailed results and discussion). At sufficiently low 
field strengths or/and high-Ko, where the magnetic torque be- 
comes weak, this analysis predicted a regime of stability due 
to magnetic diffusion. This implies that the growth timescale 
of such instability considered in this work should be signif- 
icantly lower than the dynamical timescale of the disk. The 
detailed calculation of the instability is beyond the scope of 



this work. 

New is the high-field solution found (leftmost intersection). 
By the same line of reasoning as above, this point is expected 
to be stable, since the slopes of the dashed and solid curves are 
reversed here. It is, however, somewhat outside the assump- 
tions made, since the MRJ turbulence that was assumed for 
the magnetic diffusion is probably suppressed at these field 
strengths. Instead, instability of the strong field itself is likely 
to cause its outward diffusion (as in the simulations of Stehle 
& Spruit 2001 and Igumenschev et al. 2003). To the (uncer- 
tain) extent that this process can be parametrized in terms of 
Cj/f2, the present analysis would still apply. We speculate 
that the stability of this point is actually significant, and that 
it is relevant for the experimentally observed stability of the 
strong central flux bundles in numerical simulations of accre- 
tion onto black holes. 

4. CONCLUSIONS 

We have considered the possibility that the angular momen- 
tum of an accretion disk is removed predominantly by out- 
flows driven by the accreted field (Bisnovatyi-Kogan & Ruz- 
maikin 1974, Blandford 1976). An obstacle to th is proposal 
has be en the realization (van Ballegooijen 1989, Lubo w et al.l 
(1994)) that outward diffusions of the accreted field is fast 
compared to the inward accretion velocity in a geometrically 
thin accretion disk if the value of the Prandtl number P m is 
around unity. Revisiting this problem, we find that even mod- 
erately weak fields can in fact cause sufficient angular mo- 
mentum loss via a magnetic wind to balance outward diffu- 
sion. The estimate in Eq. (f2Tb shows that, at P m = 1, a field 
with a magnetic pressure as low as ~ 10~ 3 of the gas pressure 
p c at the disk midplane has a chance of facilitating its own ac- 
cretion by driving a moderately strong magnetic outflow. This 
is due more or less to compounding numerical factors of order 
unity. In particular when MRI turbulence produces the rela- 
tively low effective viscosity a ~ 0.01 that is seen in several 
numerical simulations. 

Using a simplified model for the vertical disk structure, we 
have studied the conditions for existence of such stationary 
equilibria between wind-induced advection and outward tur- 
bulent diffusion in more quantitative detail. Two equilibrium 
points are found, one at low field strengths corresponding to 
a plasma-beta at the midplane of order several hundred, and 
one for strong accreted fields, /3 ~ 1. We surmise that the first 
is relevant for the accretion of weak, possibly external, fields 
through the outer parts of the disk, while the latter one could 
explain the tendency, observed in full 3D numerical simu- 
lations, of strong flux bundles at the centers of disk to stay 
confined in spite of strong MRI turbulence surrounding them 
(e.g. Beckwith et al. 2009). These authors also identify the 
mechanism responsible for maintenance of the bundle against 
the outward diffusion that one might expect from the turbu- 
lence surrounding it. Unlike the present model, this mecha- 
nism does not depend on the presence of a wind. 
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